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The effect of interactions on dynamics of coupled motor proteins is investigated 
theoretically. A simple stochastic discrete model, that allows to calculate explicitly 
the dynamic properties of the system, is developed. It is shown that there are two 
dynamic regimes, depending on the interaction between the particles. For strong in- 
teractions the motor proteins move as one tight cluster, while for weak interactions 
there is no correlation in the motion of the proteins, and the particle separation 
increases steadily with time. The boundary between two dynamic phases is specified 
by a critical interaction that has a non-zero value only for the coupling of the asym- 
metric motor proteins, and it depends on the temperature and the transitions rates. 
At the critical interaction there is a change in a slope for the mean velocities and a 
discontinuity in the dispersions of the motor proteins as a function of the interaction 
energy. 

Motor proteins are active enzyme molecules that are important for molecular transport, 
force generation and transfer of genetic information in biological systems They 
move along the rigid linear tracks by utilizing the energy of hydrolysis of ATP or related 
compounds, and the chemical energy is transferred into the mechanical work with a high 
efficiency. However, the mechanisms of the mechanochemical coupling in the motor proteins 
are not fully understood Q|. 

Structural and biochemical studies of the motor proteins reveal that they consist of 



many domains and subunits 



DQDD 



and frequently these subunits also have enzymatic 



activity. An example is the helicase motor protein RecBCD 



6] 



that corrects the DNA 



breaks and defects by unwinding the double-stranded DNA molecules into separate chains 
3, 0,0, Q- It has three protein subunits, of which two domains, RecB and RecD, also exist 
as independent motor proteins Experiments indicate that the complex motor protein 

RecBCD moves significantly faster than the individual RecB and RecD subunits [9j. For 
many other motor proteins the coordination between internal domains have a strong effect 



on the dynamic properties jlll . ,12j. In addition, many motor proteins work in large groups 
I2J, |3|, although the mechanism of such coordinated motion is largely unknown. In recent in 

n 

vivo experiments p^J the transport of organelles by kinesin and dynein motor proteins have 
been investigated. Although the kinesins and dyneins move in opposite directions on the 
microtubules, it was found that they do not work against each other. Apparently, the motor 
proteins moving in different directions coordinate the overall transport of the organelles. 
These experimental findings suggest that the inter-domain coupling in the motor proteins 
and the interaction between different motor proteins have a strong effect on functioning of 
these biological molecules. However, theoretical investigations of these phenomena are still 
Ihehted flQ Q. Receat.y, we — A a theoretical approach to e X p,ai, the interna! 
interactions in the motor proteins [lfij . and it was successfully applied to understand the 
dynamics of single RecBCD helicases. The purpose of this work is to investigate the general 
effect of interactions inside the motor proteins and between the molecules on the dynamic 
properties of the system. 

We assume that there are two interacting particles that move along the parallel linear 
tracks, as shown in Fig. 1. This model describes the motion of RecBCD helicases with two 
active subunits on different DNA strands |l5|, or it might correspond to the transport of 
two interacting motor proteins (kinesins, dyneins) on the parallel filaments (microtubules). 
The positions of the particles A and B are defined by integers I and m, respectively, on the 
corresponding lattices. It is assumed that the interaction between particles favor compact 
vertical configurations, while the potential energy of the non- vertical configurations is larger, 
U(l,m) = Uo + e\l — m\, where the parameter e > specifies the interactions. This potential 
of interactions seems reahstic for tee orotioa of hehca^ Q, where at each step of the leadha, 
subunit the bond between two strands of DNA should be broken, and it leads to the linear 
dependence of the interaction energy on the distance between the subunits. 

We introduce P(l, m; t) as a probability to find the system in the configuration where A is 
at the position I on the first track and B is at the position m on another track at time t. The 
dynamics of the system can be described by a set of transition rates that depend not only 
on the particle type, but also on the position of the particles. For configurations (I ± k, I) 
[k > 1], the trailing particle can move forward (backward) with a rate Uji (wji), where 
j = a or b corresponds to the particle A or B, respectively. At the same time, the leading 
particle can jump forward (backward) with a rate Uj2 {wj2}- For the vertical configurations 
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(/, I) each particle can hop forward with the rate Uj 2 or it can move backward with the rate 
Wji. see Fig. 1. Note, that in our model the transition rates do not depend on the particles 
separation k = \l — m\, but only on the "type" of transition: where it leads to a more 
compact configuration (k decreases) or a less compact (k increases). This is because of the 
linear potential of interaction, U — Uq + ek, and it leads to the energy difference between 
two consecutive configurations being equal to e, independent of the particle separation k. 
The transition rates are related via the detailed balance relations: 

^ = ^eM+c/k B T), ^ = ^eMs/k B T), (1) 

Wji Wj Wj 2 Wj 

with j = a or b, and where Uj and Wj are the hopping rates in the case of no interaction 
between the particles (e — 0). 

The dynamics of the system is governed by a set of Master equations for the probability 
distribution function P(l,m;t), 
dP(l,l;t) 



dt 



= u al P(l - 1, 1; t) + w a2 P(l + 1, /; t) + u bl P(l, l-l;t) 



+w b2 P(l, I + 1; f) - (u a2 + w al + u b2 + w bl )P(l, Z; t); (2) 

dP(l,l-k;t) 



dt 



u a2 P(l - 1, 1 - k; t) + w a2 P(l + l,l-k;t) + u bl P(l, I - 1 - k; t) 



+w bl P(l, I + 1 - k; t) - (u a2 + w a2 + u bl + w bl )P(l, l-k;t); (3) 

dP(l — k l't) 

v , ' 1 ; = u al P(l — 1 — k,l;t) + w al P(l + 1 - k, I; t) + u b2 P(l — k,l — l;t) 

dt 

+w b2 P(l - k, I + 1; t) - (« i + w a i + u b2 + w b2 )P(l - k, I; t). (4) 



-co +oo 



At all times these probabilities satisfy the normalization condition, ^2 P(l, m; t) = 1. 

l=— oo m=— oo 

The solutions of the Master equations can be found be summing over all integers / and m 
at the fixed particle separation k. Defining new functions, 

+oo +oo +oo 

Z=— oo Z=— oo i=— oo 

it can be shown then that in the stationary-state limit, 

Po,k = Po,o(Po) k , Pi,k = Po,o(Pi) k , (6) 

where 

n _U a2 + W bl _U b2 + W al 

Po — : , Pi — : • {< 

Ubl + w a2 u al + w b2 
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These auxiliary functions play a critical role in our analysis. When ft < 1 and ft < 1, using 
the conservation of probability, we obtain 

;i -#,)(! -ft- 



—(A 



0,1. 



(8) 



i - A)A 

This means that the vertical configuration (k = 0) is the most probable one, and the 
probabilities of the less compact configurations are exponentially decreasing functions of the 
particle separation k. In this dynamic phase, the particles A and B correlate their overall 
motion. From the knowledge of the stationary probabilities and the transition rates, the 
dynamic properties of the system, such as the mean velocity V and dispersion (effective 
diffusion constant) D of the center of mass, can be calculated as 

1 



V, 



CM 



PoPi 



[(Ua2 ~ PoW a2 )(l - Pi) + (U b2 ~ (3lW b2 )(l ~ ft)] 



(9) 



and 



D 



CM 



l - A)A 



1/ . o \ (A + w a2 )(u a2 - p A ) 

~{Ua2 + P0W a2 ) ■ 

2 Uh\ + w a2 

(Ai + w b2 )(u b2 - PiAi] 



+ \ -{u b2 + (3iw b2 ) - 



Ual + W b2 



;i-a)+ 

(i-A))] 



where the coefficients Ai are given by 

„ (3\{u al - U a2 ) + (3 f3iW a2 - W a i 
Aq — 



Ai 



Po(u bl -u b2 )+ p (3iw b2 - w bl 



(10) 



1 - PoPi 1 - PoPi 

The dynamic properties of the individual particles coincide with the dynamic properties of 
the center of mass of the motor protein cluster. In this case, it can be shown that the average 
distance L between the particles is always finite (in units of lattice spacings), 



PoPi 



/3 (1-A) , Pxil-Po) 



(12) 



1-Po I -Pi 

The situation is very different when, at least, one of ft > 1 (i = or 1). Then from Eq. 
© it can be concluded that less compact configurations (large k) dominate the steady-state 
dynamics of the system. In this regime the particles A and B move independently from each 
other with mean velocities (assuming A is the leading particle) 



Va = u a2 - w a2 , Vi 



B 



Ubl - w bl , 



and dispersions 



Da = (u a2 + w a2 )/2, D B = (u bl + w bl )/2. 



(13) 



(14) 
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The dynamic properties of the center of mass of the motor protein cluster is given by 

V CM = \(y A + V B ), D CM = h y D A + D B ). (15) 

Furthermore, the average particle-particle separation L is steadily increasing with time. 

The boundary between two dynamic regimes is determined by the condition Pq = 1 
and Pi < 1, or (5\ = 1 and p < 1, and it depends on the transition rates and energy of 
interaction. Using the detailed balance conditions it can be argued that the transition 
rates can be expressed as 

U ji = Ujj 1 ' 6 ' 1 , wji = Wjj"^ 1 , u j2 = Uj-i~ dj2 , w j2 = wj-f 1 ' 9 ' 2 , (16) 



where 7 = exp(e/ksT), and j = a or b. The coefficients 9ji are energy-distribution factors 
that determine the effective splitting of the interaction energy between the forward and 
backward transitions In the simplest approximation, we assume that all energy- 

distribution factors are approximately equal to each other, < 8ji « 9 < 1, because they 
describe similar transitions in the motion of the individual motor proteins jl5 1 . More general 
situation with state-dependent energy-distribution factors can also be analyzed. Substituting 
Eq. ((IH]) into the expressions (jZJ), we obtain 

Pol = {PilY 1 = K + w b )/{u b + w a ). (17) 

Then the boundary between two dynamic phases corresponds to the critical value of the 
interaction energy, 

Er = k „ T In ( ^) > 0. (18) 
\u b + w a J 

It is important to note that the critical interaction depends on temperature, and in the 
transport of the identical particles (A = B) the critical interaction is always zero. This 
indicates that the dynamic phase transition can only be observed for the coupling of the 
asymmetric motor proteins. 

The existence of two dynamic phases in the transport of interacting asymmetric motor 
proteins can be understood using the following arguments. Consider the configuration where 
the particle A is k sites ahead of the particle B and e = 0. The effective rate of the transition 
to the configurations where two particles are separated by k + 1 sites is equal to u a + w b , 
while the effective rate for k + 1 — >• k transition is given by u b + w a . The free energy change 
of making the particle configuration less compact (k — * k + 1) can be written as AG(0) = 



6 

— keTln [ Ua + Wb ) < Id, Il6|l , assuming that u a +Wb > Ub+w a . If there is interaction between 
the particles, then the free energy change increases by the value of e, AG(e) = AG(0) + e. 
The boundary between two regimes corresponds to AG(e c ) = 0, and it leads to e c = |AG(0)|. 
Thus, for strong interactions (e > e c ), it is thermodynamically unfavorable to make less 
compact configurations. The particles cannot run away from each other, and they move 
as one tightly-coupled cluster. For weak interactions (e < e c ), the favorable free energy 
change of making less compact particle configuration cannot be compensated by the energy 
of interaction. As a result, the distance between particles grows linearly with time, and they 
move in the uncorrelated fashion. 

The dynamic properties of interacting motor proteins are different in two phases, as shown 
in Figs. 2 and 3. The mean velocity of the center of mass changes the slope at the critical 
energy of interaction, while the mean velocities of the individual particles converge to a 
single value - see Fig. 2. The effect of the interaction is much stronger for the dispersions. 
As illustrated in Fig. 3, there is a jump in the mean dispersion of the center of mass at 
the phase boundary. In addition, the mean dispersions of the individual particles do not 
converge to a single value. This discontinuity in the dispersions is a clear sign of the dynamic 
phase transition in the system. 

In order to illustrate our approach, we consider a simplified model of the motion of the 
interacting motor proteins that can only step forward, i.e., w n = Wb = 0. This model seems 
reasonable for the description of RecBCD helicase transport [15j , since the experiments indi- 
cate that the backward transitions are moves faster 
than the particle B (u a > Ub), the critical interaction can be written as e c = kBTln(u a /ub)- 
For RecBCD motor proteins, where the transition rates for subunits can be approximated as 
u a = 300 and Ub = 73 nucleotides/s the critical interaction is e c ~ 1.4k B T. Theoretical 
analysis 15] estimates the energy of interaction between the subunits in RecBCD as ~ QksT, 



implying that t 



experiments 



his motor protein moves in the strong coupling regime, in agreement with 
3, sl, 0, 13 • Using Eqs. ()9I10I1 lj) . it can be shown that for large interactions 



the dynamic properties of the system are given by 

Vcm(s > e c ) = ^ + Ub ^\ D CM (e > e c ) = V CM (l - 2 ) . (19) 
1 + 7 1 V 7(1 + 7 ) J 

In the weak coupling regime, from the expressions ()13I14I15|) it can be derived that 

V CM (e < e e ) = («« + u b j)^- e , D CM (e < e c ) = Vcm/8. (20) 
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The jump in the dispersions at the critical interaction is equal to 

A/J = K/4)(V<T 

Although for this simplified model the dispersion jump is always positive, it can be shown 
in general that the discontinuity might have any sign. 

The presented theoretical analysis of the dynamics of the coupled motor proteins is based 
on the simplified picture that neglects many important features of the biological transport. 
The intermediate biochemical states, sequence dependence of the transition rates, protein 
flexibility have not been taken into account in this approach. However, it is expected that 
these phenomena will not change the main prediction of our analysis - the existence of the 
dynamic phase transitions that depend on the interaction between the particles. The most 
crucial assumption in our approach is the assumption of the linear potential of interactions. 
An important question is if the predicted dynamic phase transitions will survive for more 
realistic potentials of interaction between proteins. 

In summary, the effect of interaction between the motor proteins is investigated by an- 
alyzing explicitly a simple stochastic model. Using the explicit formulas for the dynamic 
properties, it is shown that there are two dynamic phases for asymmetric motor proteins 
depending on the interaction energy. Below the critical interaction the particles do not corre- 
late with each other, while above the critical interaction the particles move as a tight cluster. 
The origin of these phenomena is the balance between the chemical free energy change and 
the change in the energy of interactions for different transitions. The critical interaction 
depends on the transition rates and it can be modified by changing the temperature. Our 
method is applied to analyze the dynamic phase of RecBCD helicases in agreement with the 
experiments. This theoretical approach suggests a new way of investigating and controlling 
biological transport processes at the nanoscale level. 
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2< + Ub -1 



> 0. 



(21) 
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Figure Captions: 

Fig. 1. Schematic view of the motion of two interacting motor proteins. Transition rates u a % 
and w a i (i = 1 or 2) describe the motion of the particle A (small circles), while uu and wu 
are the transition rates for the particle B (large circles). Any configuration is specified by 
the integers I and m for the positions of the particle A and B, correspondingly. The energy 
of interaction in the configuration (l,m) is equal to \l — m\e > 0. 

Fig. 2. Relative velocities for the coupled motor proteins as a function of the interaction 
energy. Solid line corresponds to the relative velocity of the center of mass of the particles, 
while the dotted lines are the relative velocities of the individual particles below the critical 
interaction. The parameters used for calculations are: u a — 4, w a — 0.1, Uf, — 1, iu& = 0.1 
and 6 = 0.02. 

Fig. 3. Relative dispersions for coupled motor proteins as a function of the interaction 
energy. Solid line corresponds to the relative dispersion of the center of mass of the particles, 
while the dotted lines are the relative dispersions of the individual particles below the critical 
interaction. The parameters used for calculations are the same as in Fig. 2. 
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Figure 1: E. Stukalin, A. Kolomeisky, Physical Review Letters. 
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Figure 2: E.B. Stukalin, A.B. Kolomeisky, Physical Review Letters. 
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Figure 3: E. Stukalin, A. Kolomeisky, Physical Review Letters. 



